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Half a century ago Brian Josephson made 
a series of striking predictions related to 
a tunnelling barrier sandwiched between two 
superconductors.^ One particular prediction, 
later became known as the a.c.-Josephson effect, 
said that under a finite d.c. bias the tunnelling 
current contains an a.c. supercurrent component, 
oscillating at microwave frequency. This predic- 
tion was experimentally verified through obser- 
vation of the junction current- voltage character- 
istics modified by the interaction with the elec- 
tromagnetic radiation, and direct detection of 
the microwave radiation emitted from the junc- 
tion itself.^ It had also been established that 
the behaviour of the d.c. current-voltage char- 
acteristic is determined by the high-frequency 
dynamics of the Josephson junction,^ and can 
be used to study various systems such as atoms 
in the electron-spin resonance,^ optical phonons 
in high-Tc superconductorsf^ and vibrating Nb 
molecules.^ Here we present an InAs nanowire 
Josephson junction device^iiS where a vibrating 
nanowire weak link plays the role of a nanoelec- 
tromechanical resonator. The fiow of the a.c.- 
Josephson current through the junction enabled 
actuation and detection of vibrational modes of 
the resonator by means of simple d.c. transport 
measurements. 

There are two regimes a Josephson junction can be 
operated ini The first, the so-called d.c. -Josephson ef- 
fect, is characterised by a d.c. superconducting current 
Is = Ic sin (A^), where Ic is the critical current, and Acf) 
is the superconducting wave functions phase difference. 
In the second regime corresponding to the a.c.-Josephson 
effect the tunnelling current exceeds the value of /c, and 
consists of two components. Apart from a normal dissi- 
pative current manifested by the voltage drop the to- 
tal current contains the oscillating superconductive part 
Is = /c sin (cjjt -h ^o), where 0o is the initial phase dif- 
ference. In this case the phase difference increases lin- 
early in time t with the rate set by the d.c. voltage as 
cjj = {2e/h)V. 

Because of the high frequency of the a.c.-Josephson 
current {2e/h ^ 483.6 MHz//iV), the shape of the 
Josephson junction current- volt age (/ — V) characteris- 
tics is determined by its electrical impedancei^i^ One of 
the most vivid demonstrations of this fact was the ex- 
periment of Coon and Fiske,^ where the I — V character- 
istics of the Josephson junction coupled to the electro- 
magnetic resonator, exhibited constant- volt age step fea- 
tures (self- induced Shapiro steps). In this experiment 
the steps occurred when the voltage V was such that the 



a.c.-Josephson oscillations were in resonance, namely, at 
voltage Vn = hujn/2e, where uOn is the frequency of the n- 
th resonator mode. This work illustrated how a Joseph- 
son junction coupled to an external resonant system can 
be simultaneously used as a source and detector of the 
microwave radiation. This remarkable property allowed 
the use of a Josephson junction as an ultra-small power 
spectrometer to study atomic and vibrational spectraA^ 
Recently, the combination of the Josephson junction 
and nanoelectromechanical system^^ (NEMS) has at- 
tracted increasing attention. It was proposed to use 
the nanomechanical resonator to store the Josephson 
junction phase qubit states ^li^ and to build a co- 
herent source of microwave radiation based on an ar- 
ray of Josephson junctions^^ In addition, a vibration- 
modulated Josephson junction has been predicted to 
exhibit self-induced Shapiro steps in the I — V char- 
acteristics at a frequency of mechanical resonance^ 
and the weak-link made of a doubly-clamped nanowire 
can be used to pump the supercurrent into nanowire 
vibrations. In all these aspects NEMS act similarly to 
an external microwave resonator coupled to a Josephson 
junction. 

Here we report observation of the a.c.-Josephson oscil- 
lations being in resonance with the longitudinal acoustic 
wave of an InAs nanowire. Using a previously devel- 
oped technique^- we fabricated a new type of nanowire 
Josephson junction device^iiSl (Figs. [1^ and [TJd). In 
our field-effect type device the 240 nm-long nanowire 
weak-link is suspended in vacuum and clamped by two 
superconducting Ti/Al contacts, making the nanowire- 
based NEMS^^ an integral part of the Josephson junc- 
tion. The transport experiments with the current-driven 
Josephson junction described here were carried out in- 
side the dilution refrigerator with Tbase ^ 15 mK. All 
measurements were taken at a fixed voltage on the local 
and global backgates. Figure [It shows the experimen- 
tal / — V curve typical of an underdumped Josephson 
junction; a non-dissipative branch limited by the critical 
current Ic ^ 6.5 nA, some hysteresis at around /c, and 
the resistive part, with the well-pronounced step-like fea- 
tures. To look into the details of the I — V curve we 
measured the differential resistance dV/dl as a function 
of applied voltage V (Fig. [T]i). At the subgap voltages 
(e|V| < 2 A, here 2 A is the superconducting gap in the 
Al contacts) the differential resistance reveals an over- 
all decrease with some dips originating from the excess 
current of the so-called Multiple Andreev Reflection^ 
(MAR). The dips in the resistance occur when the super- 
conducting gap 2A is an integer multiple of the applied 
bias, which defines the position of the MAR signature 
as Vm = =b2A/em, where m is integer. In our experi- 
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FIG. 1. The suspended InAs nanowire-based Josephson junction and its / — characteristics, (a) A schematic 
representation of the suspended nanowire device and the experimental set-up. The InAs nanowire (NW) is suspended over 
a trench predefined in Si02, and clamped by two Ti/Al contacts. An in-plane magnetic field B was applied during the 
experiment, and it was only few degrees off the nanowire axis, (b) The top view of a typical device obtained by a scanning 
electron microscope. The scale bar is 200 nm. The fake colours indicate different parts of the device. Blue is the Ti/Al 
contacts, yellow is the nanowire, and gray is the local backgate below the nanowire. (c) The I — V characteristics of the 
nanowire Josephson junction taken at = 12.3 V. (d) The differential resistance dV/dl as a function of the voltage V taken 
at the same Vg as in (c). The red arrows show first four orders of MAR at voltage Vm- The position of the first MAR signature 
(m = 1) gives the superconducting gap 2A ^ 230 /xeV and the estimate for Tc = 2A/(3.52A;b) ~ 760 mK. The black arrows 
mark the resistance resonances caused by the acoustic wave actuated by Josephson oscillations at voltage Vn. 



ment we clearly observed MAR signatures with m up to 
four (red arrows in Fig. [T]i). However, at lower biases, 
the subgap features become sharp well-recognised asym- 
metric resistance resonances (black arrows in Fig. [T]i), 
corresponding to the steps on the I — V curve (Fig. [T]^). 
The distinct nature of these resonances becomes appar- 
ent from the gray-scale plot of dV/dl measured at dif- 
ferent magnetic fields shown in Fig. [2^. The position of 
the MAR signatures is tightly bound to the value of the 
magnetic field- and temperature-dependent gap 2 A, and 
shifts toward V = as the magnetic field 5, or tem- 
perature T, approach their critical values. As seen from 
Fig. [2^, the position of the MAR signatures for m = 1, 2 
follow the prediction of the Landau- Ginzburg theory^^ 
(red dashed curves). Whereas, the position of the re- 
sistance resonances (black arrows) remains unchanged, 
and the resonances simply vanish for B ^ Be. Quali- 
tatively the same behaviour is observed for dV/dl mea- 
sured at different temperatures. Moreover, the position 
of the resonances is independent of the backgate voltage 
(Fig.lSt), which rules out the superconducting bound 
states21i2^ as possible explanation for there origin. 

Such a peculiar behaviour indicates that these features 
are due to the self-induced Shapiro steps of the junction 
coupled to an external resonator The first resonance 
appears at around V ^ 20 /iV, which corresponds to an 



a. c. -Josephson frequency of about 10 GHz. In our exper- 
iment there was no deliberately designed microwave res- 
onator coupled to the nanowire junction. The possibility 
of a microwave cavity unintentionally formed inside our 
cryogenic setup was excluded by reproducing the same 
result in a different cryogenic system. The distributed 
parameter resonators, such as coplanar striplines, are 
ruled out by the fact that the dimensions of the de- 
vice including the leads do not exceed 1 mm, and it is 
much smaller than the radiation wavelength (~ 3 cm). 
We also considered a lumped LC circuit, where the ca- 
pacitance is defined by the area of the device bonding 
pads (~ 3 pF), and the inductance is defined by the de- 
vice leads (~ 1 nH). The estimated resonant frequency 
(~ 3 GHz) appears to be well below the measured one. 

The important clue to the understanding of the ob- 
served resistance resonances lies in the fact that the 
nanowire junction is suspended and its possible mechan- 
ical oscillations are not dumped by the substrate. In- 
deed, a quick estimation shows that the a. c. -Josephson 
frequency of the resonances turn out to be close to the 
frequency of the first two longitudinal acoustic modes of 
the doubly-clamped InAs beam^i /n = {'d^/^L) xn^ 9A 
GHzxn, where n - is integer, i^g, = 4.5 x 10^ m/s is 
the speed of sound in wurtzite InAs nanowire along the 
(111) direction, and L = 240 nm is the length of the 
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FIG. 2. The magnetic field- and gate voltage-independent position of the resistance resonances, (a) The gray-scale 
plot of dV/dl made in the V — B plane at T = Tbase- The red dashed curves show the position of the first two orders of MAR 
expected from the Landau-Ginzburg theory,— A(5)/A(0) — y^l — {B/BcY^ where Be = 63 mT is the critical field. The black 
arrows mark the magnetic field-independent position of the resistance resonances. The position of the resistance resonances 
stays unchanged while the peak amplitude decreases while B approaches its critical values Be. (b) The gray-scale plot of 
dV/dl made in the V — Vg plane at the same temperature as in (a). The red arrows point at the first four MAR signatures 
seen in Fig. [TJl. The black arrows show the gate voltage-independent resistance resonances. The horizontal red dashed line 
marks the dV/dl trace shown in Fig. [TJl taken at Vg = 12.3 V. The changing shades of gray quasi-periodic in Vg are due to the 
Fabry-Perot quasiparticle conductance oscillations.— 



beam. This suggests that the resistance resonances are 
caused by interaction of the a.c.-Josephson current with 
the acoustic modes of the nanowire. Since the mode fre- 
quency is defined only by the nanowire length the posi- 
tion of the resonances caused by them should be inde- 
pendent of the values of A and Vg^ which is confirmed by 
our observations. 

To further understand this phenomenon we used a sim- 
ple model of the Josephson junction coupled to an InAs 
nanowire-based Bulk Acoustic Resonator— (BAR). 

The junction is described by the RSJ-model^ (Fig.[3ti), 
where an ideal junction J is connected in parallel to the 
normal state resistance i?n, and the total junction capac- 
itance Cy As mentioned before, the shape of the I — V 
characteristics is governed by the high-frequency prop- 
erties of the junction, namely by the relative impedance 
of Rn and Cj, and characterised by the parameter /3 = 
cJci^nCjr where Uc = {2e/h)IcRn is the lowest frequency 
of the Josephson oscillations. If Cj is large {f3 » 1) the 
junction is 'underdumped' and the I — V curve is linear. 
In the opposite limit of small Cj (/3 << 1), the junction 
is 'overdumped' and the I — V curve is parabolic. The 
experimental values of the critical current (/c = 6.5 nA) 
and the normal-state resistance {Rn ~ 2 kOhm) gave us 
the lowest value of the Josephson oscillations frequency 
a;c/27r ^ 6.3 GHz, and /3 ^ 30 (with Cj ~ 3 pF), indicat- 
ing that the junction is 'underdumped'. 

In our device capacitance Cj is a result of the parallel 



connection of the leads capacitance and the piezoelec- 
tric capacitor formed by the wurtzite InAs nanowire^ 
In the case of a piezoelectric capacitor the external a.c. 
electric field excites acoustic waves inside the dielectric 
media and vice versa. The process of energy conver- 
sion from the electrical to the mechanical domain pro- 
foundly affects the capacitor impedance, and is the basis 
of the BAR operation. According to the one-dimensional 
Mason model^^ the total effective impedance Zq{uj) of a 
piezoelectric capacitor is 



1 



tan 



ujL 



(1) 



where uo is the angular frequency, is the piezoelectric 
coupling constant, is the speed of sound in the piezo- 
electric media, and L is the thickness of the piezoelectric 
material. It is clear from Eq. ([1]) that the longitudinal 
acoustic standing wave with uo^ = {7r/L)i9a causes the 
total impedance to diverge {Zq{uJj.) oo). 

Now we assume that instead of Cj the Josephson junc- 
tion dynamics is defined by the impedance Zq of the InAs 
nanowire BAR (Fig. [3)3). While the BAR is off resonance, 
the I — V curve follows the one expected for the under- 
dumped regime (blue curve. Fig. [Sh). When the volt- 
age V is such that the Josephson oscillations are in reso- 
nance with the BAR its impedance increases, forcing the 
Josephson junction to cross from the underdumped to the 
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FIG. 3. The RSJ-model and I — V characteristics of the Josephson junction coupled to the Bulk Acoustic 
Resonator, (a) The RSJ-model circuit diagram of a typical Josephson junction, (b) The RSJ-model of the Josephson 
junction coupled to the BAR. (c) The schematic plot of the Josephson junction I — V characteristics. If the capacitance Cj is 
large (/3 ^ 1) the a.c- Josephson current is shunted by Cj, and its time-averaged contribution to the I — V characteristics is zero. 
Thus for I > Ic the I — V characteristics is of a normal conductor V — IRn (black dashed line). In the opposite case of small 
Cj (/3 <^ 1) the alternating current is forced to flow through Rn causing the a.c. -Josephson current to be frequency- modulated. 
The time- average of th e frequency-modulated a.c. current is non-zero, which makes the total I — V characteristics parabolic^ 
V/Vc = V(///c)2 - 1, where K = Rnh (blue curve). The I — V characteristics of the underdumped Josephson junction 
coupled to the BAR is shown by the red curve. When the voltage V = Vr the Josephson oscillations actuate the acoustic 
resonance in BAR (Ze(cj) oo), which forces the junction to cross from the underdumped to the overdumped regime, and 
results in the current peak SI. For simplicity we assumed that the resonance in Ze has a Lorentzian lineshape. In the case 
of the current- driven junction, the current peak takes a step-like shape with a multi-valued region showed by the red dotted 
curve. The green curve on the right-hand side is the corresponding dV/dl as a function of V. (Inset) The actual I — V curve 
measured at T = Tbase and B = 0. (d) The zoomed-in plot of the dV/dl resonance taken at different temperatures and zero 
magnetic field. Identically to the magnetic field dependence (Fig. [2^), the position of the resonances is independent of the 
temperature, which rules out their MAR-related origin. We used the interception point of curves taken at different T to identify 
the resonant voltage V^n- The vertical dashed lines designate the voltage Vn which corresponds to the first three longitudinal 
harmonics expected for the doubly-clamped InAs beam with length L = 230 nm and speed of sound i9a — 4.5 • 10^ m/s. The 
respective frequencies are given in brackets. 



overdumped regime (red curve, Fig.[3t). This crossing re- 
sults in a cusp-like increase of the d.c. -current through 
the junction 51 with the maximum at K = hu^jle^ 
where uj^ is the mechanical resonance frequency. In the 
case of the current-driven Josephson junction the cur- 
rent peak on the I — V curve becomes a step-like feature 
and corresponds to the asymmetric resonance in dV/dl 
(green curve, Fig. [3b). 

Using this model, we determined the resonant fre- 
quency of the nanowire BAR actuated by the Joseph- 
son oscillations. Two resistance peaks shown in Fig. [3]i 
correspond to the resonant coupling of the first two lon- 
gitudinal acoustic modes of the doubly-clamped InAs 
nanowire, with onset at the voltage determined as 

" ~ 2e ~ 2e ' 2L* ' 

where n = 1, 2 and L* = 230 nm is the length of the 
acoustic resonator (i^a = 4.5 x 10^ m/s). Note, the length 
of the resonator L* found here is close to the physical 
length of the nanowire Josephson junction (240 nm). De- 
spite the fact that the resonance in Zq requires a non-zero 
spatially averaged strain^ which forbids the appearance 



of even harmonics, we have observed the second harmonic 
as well. The reason for that is the slightly asymmetric 
geometry of our device. 

The piezoelectric coupling of the Josephson oscillations 
to the acoustic mode of a crystal resonator has been 
theoretically studied for the case of an array of Joseph- 
son j unctions From this point of view our experimen- 
tal realisation can be seen as a special case of a single 
weakly coupled junction. Other possible coupling mech- 
anisms such as the magnetomotive forced and modu- 
lation of the junction dimension^ are ruled out since 
the resonance is observed at 5 = 0, and the junction 
length L is fixed. It is worth mentioning that the strong 
coupling between the longitudinal vibrations and single 
charge tunneling by means of electron- vibron coupling2S 
has been observed in suspended carbon nanotubes j^^i^^ 
However, the underlying physical principals responsible 
for the electron- vibron coupling is entirely different from 
the one studied in this work. 

In conclusion, we fabricated an InAs nanowire nano- 
electromechanical Josephson junction, and observed the 
resonant coupling of Josephson oscillations to the lon- 
gitudinal acoustic modes of the nanowire. The reso- 
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nant coupling reveals itself as voltage steps in the junc- 
tion I — V characteristics occurring when the Joseph- 
son oscillations are in resonance with the acoustic modes. 
We showed that the Josephson junction coupled to the 
NEMS can be used as a self-actuating probe to study 
vibrational states. The advantage of this tool is in 
its simplicity, the NEMS is probed by routine d.c- 
measurements, with no need for external high-frequency 
equipment Jiii^ii^ Also, our experiment demonstrates the 
possibility of building a nano-scale source of coherent mi- 
crowave radiation based on an array of Josephson junc- 
tions coupled to a high-Q piezoelectric resonator 

METHODS - The high-quality InAs nanowires used 
in this work have been synthesised by the Molecular- 
Beam Epitaxy in Riber 32 solid source system by the 
Au-assisted Vapour-Liquid- Solid (VLS) method. The 
nanowires were grown along (111) direction on a (011)- 
oriented InAs substrate at 400 °C and the group V/III 
ratio of 100. All nanowires have a pure wurtzite crys- 
tal structure with only one stacking fault per several 
nanowires 

The basis of the suspended nanowire Josephson junc- 
tion was a p+— Si wafer with 300 nm thick Si02 ther- 
mally grown on top. First, an array of 250 nm x 20 /am 
trenches was chemically etched in the Si02 to depth of 
about 150 nm. This was followed by evaporation of the 
local metal backgate on the bottom of the trenches. Then 
nanowires were randomly distributed across the prede- 
fined trenches so that a small segment of a nanowire 
is suspended over the trench. E-beam lithography was 
employed to pattern Ti/Al contacts (5 nm/100 nm) 
on top of the parts of the nanowire supported by the 
banks of the trench, leaving the suspended segment un- 



touched (Figd^). To ensure high contact transparency 
the nanowire native oxide was removed with the 0.2% 
solution of (NH4)2Sa: {x = 1.5 M) prior to evaporation of 
the contacts. The nanowire junction was about 240 nm 
long with a diameter of about 50 nm (Fig[T]3). The ex- 
perimental data obtained from another device is given in 
the Supplementary Information. 

All transport measurements were carried out inside a 
^He -^He dilution refrigerator with heavily filtered signal 
lines and Tbase ^15 mK. The I — V characteristics and 
differential resistance as a function of d.c. voltage were 
measured in a quasi four-terminal constant current con- 
figuration. The voltage drop across the device was ampli- 
fied with a low-noise voltage preamplifier and registered 
using a high-speed digital voltmeter or lock-in amplifier 
(Fig. [TJi). The a.c. excitation current /ac was set to be 
below 400 pArms, so that the excitation voltage drop 
across the sample was not greater than /cbT. 
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I. CONTACT TRANSPARENCY 
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FIG. 1. The excess current of the In As nanowire Josephson junction described in the main text taken at Vg = 12.3 V. 



To assess the transparency of the superconducting Al contacts made to the InAs nanowire we measured the so-called 
excess current^ /exc- The excess current of the Josephson junction is related to the presence of Andreev Reflections, 
and defined by the average strength Z of the interface barrier between the normal weak-link and the superconducting 
contact. The value of /gxc can be determined as an intersection of the I — V curve taken for V > ±2A/e, and 
extrapolated down to V = as shown in Fig. [H Using this method we found that /exc ^ 38 nA at = 12.3 V and 
the normalised excess current ^^^""^ to be about 1.15, which corresponds to Z ^ 0.3 (see Fig. 6 in Ref. [l|). Hence, 
the contact transparency defined as (1 + is about 90%. 



II. MAGNETIC AND TEMPERATURE DEPENDENCE OF THE CRITICAL CURRENT 




FIG. 2. Additional data on the magnetic field and temperature dependence of dV/dl and h for the nanowire device described 
in the main text. 



Figured^ shows the zoomed in plot of dV/dl as a function of V at different magnetic field shown in Fig. 2a of 
the main text. The background resistance increasing with the magnetic field made it difficult to reliably extract the 
position of the resonances. However, it is clear from the plot that the resonance position is virtually independent of 
the magnetic field. Figure [2]3 is the plot of the critical current /c as a function of the magnetic field B measured at 
T = Tbase and Vg = 12.3 V, and Fig. ^ is the critical current Ic as a function of the temperature T measured at 
5 = and K = 12.3 V. 
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FIG. 3. Simple model of the underdumped Josephson junction coupled to an external electromagnetic 
resonator^, (a) The theoretical I — V curve (blue) and experimental data (red circles) taken at T = Tbase- (b) The 
theoretical differential resistance dV/dl as a function of V derived from the junction I — V characteristics simulated for differ- 
ent values of Ic- The plot clearly shows that the position of the resonance marked by vertical dashed lines corresponds to the 
intersection point of curved taken at different Ic- 

III. SIMPLE MODEL OF THE JOSEPHSON JUNCTION COUPLED TO AN EXTERNAL RESONATOR 

To verify the assumption made in the main text that the position of the resonance is given by the interception 
of the curves measured at different temperature we used a simple theoretical model of the underdumped Josephson 
junction weakly coupled to an external electromagnetic resonator^ Here the theoretical I — V curve (Fig. [3^) is the 
superposition of the I — V curve of the under dupmed Josephson junction V = Rnl and three current peaks SI with 
Lorenzian lineshape 6I{Vr) = 0.5z/c(l + ^^)~^, where z ~ 0.5 is the coupling coefficient, ^ = 2Q{V/Vr — 1), and Q ^ 3 
is the quality factor. These peaks are caused by the resonant coupling of the a. c. -Josephson oscillations to the external 
resonator characterised by the resonant voltage and the quality factor Q. To account the effect of the temperature 
and magnetic field we assumed that the only way these parameters affect the I — V curve is through changing critical 
current Ic as illustrated in Fig. [JJd and Fig. [2t. Figure [Sb shows the differential resistance dV/dl derived from the 
theoretical I — V curve. It is clear that the position of the resonance Vr corresponds to the interception point of the 
curves with different /c, which confirms our assumption. 

It is important to note that the Bulk Acoustic Resonator has two two distinct resonant frequencies.^ The first one 
is at C(Jr = TT'i^a/^^ whcrc Z(co'r) oo, which corresponds to the resonance of a parallel LC-circuit. Another one occurs 
at cjr = 'i^a/^ [tt^ — S/c^] , whcrc Z(cJr) 0, and it corresponds to the resonance of a series LC-circuit. The later one 
was not observed in our experiment, presumably, because the Josephson junction is already strongly underdumped 
by the relatively large value of Cj, and the additional shunting of the a. c- Josephson current does not take any visible 
effect on the I — V curve. 



IV. SPATIALLY AVERAGED STRAIN S IN THE CASE OF AN ASYMMETRIC NANOWIRE DEVICE 

As mentioned in the main text, the non-zero spatially averaged mechanical strain forbids the appearance of the 
even resonant modes of BAR. Indeed, the displacement current through the piezoelectric capacitor is /d = — 
(ess/e).?],— where 633 is the piezoelectric constant, e is the dielectric constant, S = {dU{z)/dz) is the spatially 
averaged mechanical strain {U{z) is displacement, z is the coordinate along the piezoelectric media), and the dot 
denotes the time derivative. The second term in the expression is the piezoelectric contribution, which is non-zero 
when the media accommodates an odd standing acoustic wave. Assuming that the displacement U{z) oc cos(/cn2;)^, 
where = ttti/L, its spatial average is defined as 




thus this integral is non-zero only for odd standing waves (odd values of n). 

However, due to the self aligned geometry of our suspended nanowire device the distance / between ohmic contacts 
is smaller than the distance L between clumped nanowire ends and the position of their centres is slightly shifted 
with respect to each other (Fig. 3^). This shift (~25 nm) results from the finite accuracy of the e-beam lithography 
alignment. In this situation the averaging with integral ([1]) should take place only between the ohmic contacts along 
the length / (from zi to Z2) as shown in Fig. IDd. For purely illustrative purpose we assume that dU/dz{z = 0, L) = 0. 
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FIG. 4. The asymmetric geometry of the experimental nanowire device and the non-zero spatiahy averaged strain S. (a) The 
scanning electron microscope image of a typical device, which shows the difference in mechanical L and electrical / length of 
the nanowire device. Also, the asymmetry due to the lithography misalignment is clearly visible. The scale bar corresponds to 
200 nm. (b) The schematically represented device cross section and the mechanical strain distribution along the nanowire for 
the first two acoustic modes in the case of the symmetrically positioned contacts. The shaded areas under the strain distribution 
curves give the value of the integral ([1]). (c) The same as (b), but with the contacts being asymmetrically positioned with 
respect to the centre of the nanowire (shifted by A^;). The resulting non-zero values of the integral ([1]) is shown for the case of 
n=2. 



For the symmetric arrangement of the contacts the first mode gives the non-zero value of S and zero value for the 
second mode. In the asymmetric case (Fig. Ht), the first mode still gives the non-zero S. Apart from that, the 
difference Az between the position of the actual nanowire centre and the centre of the contacted segment makes the 
spatially averaged strain for the second mode non-zero as well. A similar approach is used to activate even acoustic 
resonances in the so-called Composite Bulk Acoustic Resonators — 





FIG. 5. The scanning electron microscope images of the device described in the main text (a), and additional device described 
here (b). The scale bar corresponds to 500 nm. During the preparation for the SEM imaging both of the devices were broken 
down by an unintentional electrostatic discharge. The missing segments of the nanowire are marked by the red lines. 



V. ANOTHER DEVICE USED IN THE EXPERIMENT 



Apart from the device described in the main text and shown in Fig. [5^ we also measured a more complex nanowire 
configuration, namely two thinner (30 nm in diameter) InAs nanowires held to each other side by side by the Van der 
Waals forces as shown in Fig. [5]3. Both of the nanowires were suspended across the predefined trench and contacted 
in parallel with Ti/Al contacts, so they form a single weak link of the Josephson junction. 
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FIG. 6. Additional experimental data obtained from a double beam nanowire device shown in Fig. [5]3. (a) The differential 
resistance dV/dl as a function of the voltage V taken at Vg = V and T = Tbase- The red arrows show first three orders of MAR 
at voltage Vm. The position of the first MAR signature {m = 1) gives the superconducting gap A ^ 115 /xeV and the estimate 
for Tc = 2A/(3.52 /cb) ~ 760 mK. The black arrows mark the resistance resonances caused by the acoustic wave resonance in 
double beam system, (b) The gray-scale plot of dV/dl made in the V — B plane at T = Tbase- The red dashed curves show 
the position of the first two orders of MAR expected from the Landau-Ginzburg theory,^ A(B)/A(0) = y^l — (B/Bc)^, where 
Be = 63 mT is the critical field. The black arrows mark the magnetic field-independent position of the resistance resonances. 
The position of the resistance resonances stays unchanged while the peak amplitude decreases while B approaches its critical 
values Be- (c) The zoomed-in plot of the dV/dl resonance taken at different temperatures and zero magnetic field. Identically 
to the magnetic field dependence shown in (b), the position of the resonances is independent of the temperature. The vertical 
dashed lines mark the voltage which corresponds to the first three longitudinal harmonics expected for the double InAs beam 
system with length 230 nm, diameter 30 nm and speed of sound i^a — 4.5 • 10^ m/s. The coloured scale bar refers only to the 
dV/dl traces taken at different temperature. (Inset) The schematically shown submodes of double beam system. The colour 
scale corresponds to the value of the spatially distributed strain inside each of the beam (red is maximum, blue is minimum). 
It is clear that the peak in the resistance appears to be closer to the resonance of the submode given by the antisymmetric 
combination of the main mode. 



The differential resistance dV/dl as a function of the bias V taken at = V and T = Tbase is shown in Fig. [6^1. 
As can be seen from the plot, apart from the MAR signatures with m = 1, 2, 3 (red arrows in Fig. [6^) the dV/dl 
data also shows three closely spaced peak-like features in the range of V between 30 /iV and 70 /iV (black arrows in 
Fig. [6^). Similarly to the single nanowire case described in the main text, the position of these features is independent 
of magnetic field and temperature as shown in Fig. [6)3 and[6t. To evaluate the frequency of mechanical resonances of 
such double beam system we used numerical simulation made with COMSOL Multiphysics Package. We found that 
every acoustic mode turns out to be splitted into two submodes each of which is a symmetric or antisymmetric linear 
combination of a single beam acoustic mode (inset in Fig. [6t). Comparing the obtained frequencies with the position 
of the peaks in dV/dl we concluded that their position corresponds to the antisymmetric submodes for n = 1, 2, 3. 
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However, a more detailed theoretical study is needed to understand why only antisymmetric combination of the main 
modes are visible in the experiment. 
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